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Abstract. The major multi-epoch VLB A programs are described and dis- 
cussed in terms of relativistic beaming models. Broadly speaking the observed 
kinematics are consistent with models having a parent population which is only 
mildly relativistic but with Lorentz factors extending up to about 30. While the 
collimation and acceleration appears to mainly occur close to the central engine, 
there is evidence of accelerations up to 1 kpc downstream. Generally the motion 
appears to be linear, but in some sources the motion follows a curved trajectory. 
In other sources, successive features appear to be ejected in different directions 
possibly the result of a precessing nozzle. 

The launch of GLAST in 2008 will offer new opportunities to study the re- 
lation between radio and gamma-ray activity, and possibly to locate the source 
of the gamma-ray emission. VSOP-2 will give enhanced resolution and will 
facilitate the study of the two-dimensional structure of relativistic jets, while 
RadioAstron will provide unprecedented resolution to study the fine scale struc- 
ture of the jet base. 



1. Introduction 

As discussed earlier in this conference by Dave Meier, the theoretical framework 
for the formation of relativistic jets due to the accretion of matter onto a su- 
permasive black hole is well developed. Observations of jet kinematics obtained 
by multi-epoch high resolution imaging provides tests of theoretical models and 
can constrain the parameters of the jet formation. However, the interpretation 
is not straightforward. The observations do not give instantaneous multi-epoch 
snapshots of the source structure. Rather, due to relativistic effects, different 
parts of the source are observed at different times, giving the illusion of super- 
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luminal motion, the distortion of angles and misrepresentation of luminosity. In 
particular, the observed speeds and luminosities of jets are highly biased due 
to Doppler boosting within a narrow cone, and they are not representative of 
the intrinsic distributions. Finally, we note that apparent speeds derived from 
multi-epoch VLBI observations refer to the observed motion of recognizable fea- 
tures, typically referred to as "components," which do not necessarily reflect 
the underlying bulk relativistic flow, but which instead may simply reflect the 
propagation of forward and reverse shocks within the jet flow. The observations 
are further complicated since these features do not remain stable with time, but 
may break up into s ub component s , whil e other features may merge into a single 
complex feature e.g, iKadler et al.l ( 20081 ). 



2. The major multi-epoch programs 



Prior to the completion of the VLBA, observations of jet kinematics were based 
on infrequent multi-epoch VLBI observations using ad-hoc VLBI arrays gen- 
erally with too few antennas to obtain the high-sensitivity and high dynamic 
range imaging needed to unambiguously trace jet motions. A number of VLBA 
programs have now been underway for up to a decade or more at wavelengths 
between 7 mm and 13 cm and are providing new insight into jet kinematics. 
This paper reports on global source kinematics primarily obtained from 2 cm 
VLBA observations, but we also discuss the multi- wavelength kinematics of a 
few particularly interest ing sources. 



Britzen et al. ( 20071 ) have combined their 6 cm multi-epoch observations of 



compact radio sources with pre- VLBA observations to study the kinematics of 
266 radio source s . Typ ically 2 to 4 epochs were observed for each source. 

iPiner et al. (2007) have used the Radio Reference Frame Image Database 
(RRFID) to determine source kinematics for 87 sources at 4 cm over the period 
1994 to 1998. This is part of a database originally obtained for precise astrometry 
and to study earth motions and includes multi-epoch observations at both 4 and 
13 c m for up to 500 sources. 



Kellermann et al.l (|2004! ) and lCohen et al.l (|2007l ) have reported 2 cm VLBA 
studies for 110 AGN jets covering the period 1994 to 2001. Ros et al. and Lister 
et al. (in preparation) extend the observations to 2003 and 2007 respectively. 
Since 2002, the MOJAVE program (Monitoring of Jets in Active Galactic Nu- 
clei) has further extended this program to include the kinemat ics of a complete 
samp le of the the brightest 133 sources in the northern sky (jLister &: Homan 
2005) as well as 59 other sources of particular interest. Since 1994, these pro- 
grams have observed as many as 30 epochs for some sources, permitting for the 
first time clear evidence for non radial motions, bend s, accelerations and dec eler- 
ations. The MOJAVE program also includes linear ( Lister Sz Homan 20051 ) and 
circular ( Homan Sz Lister 20061 ) polarization measurements which can be used 
to investigate three-dimensional magnetic structure. The MOJAVE program is 
a large VLBA program and images of new and past 2 cm VLBA observations 
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are made available on the MOJAVE web sit^H within a few weeks following the 
correlation of the data. 

In a series of papers, the Boston University group has reported on 0.7 and 
1.3 cm VLBA observations (iJorstad et al.1 1200a . |2007| )FI These observations 
were made more frequently, typically monthly or bi-monthly, and consequently 
include a much smaller number of sources than observed by the longer wave- 
length programs. Monthly observations of 29 sources at 7 mm are continuing 
with special attention to radio counterparts of EGRET gamma-ray sources. 

Recently TAN AMI, (Tracking Active Galactic Nuclei with Australian-South- 
African Milliarcsecond Interferometry) a new VLBI program to monitor extra- 
galactic jets in a sample of sources below -30 degrees declination began to observe 
in November 2007 at 8.4 GHz and 22 G Hz with the Australian South- African 
Long Baseline Array ( Kadler et al. 20081 ). It is planned to expand this sample 



to include GLAST detected sources. 



3. Jet Properties 



Due to relativistic beaming the observed brightness temperature is enhanced by 
the Doppler factor (T Q b s = 5 x Tint). The maximum brightness temperature or 
lower limit that can be determined by VLBI depends only on flux density, mea- 
surement accuracy and baseline length, and is independent of the wavelength 
used. See Kovalev (these proceedings). For a source of ~ 1 Jy and eart h base - 
lines, this limit is ~ 10 13 K. Using the VLBA at 2 cm, iKovalev et aD <|2005T> 
have o bserved brightness temperatures up to ~ 5 x 10 13 K, while iHoriuchi et al. 
( 2004! ) have reported values of Tj, up to ~ 10 14 K using HALCA at 6 cm. This 



suggests Doppler factors up to a few hundred if T; nt is limited by Inverse Comp- 
ton cooling to value s ~ 10 n ~ 12 K. How ever, the highest observed Lorentz factors 
are only 35 to 40 ( Cohen et al. 20071 ) . much less than required to explain the 
highest observed brightness temperatures, unless these high observed values of 
T t, s are only tran s ient. 

Homan et al. ( 20061 ) have compared the distribution of observed brightness 



temperatures with relativistic beaming models and find that in their medium 
to low brightness state, most sources have T- m t ~ 3 x 10 10 K or close to an 
equipartition value. However, when they are in a flaring state, T int > 2 x 10 11 K, 
which is close to the IC limit. This correponds to an excess of particle over 

mag netic density of as much as 10 5 or more. 

Kellermann et all (2004) and IJorstad et al.l ( 20051 ) have used the charac- 



teristic variability time-scales to estimate intrinsic source size and brightness 
temperature, and find that the distribution of observed speeds and brightness 
tem perature are c onsist ent with Tj n t being close to the equilibrium value. 



Cohen et al. ( 20071 ) have shown that the distribution of observed speed and 



luminosity is consistent with relativistic beaming models having a maximum 
Lorentz factor ~ 32, and intrinsic luminosity ~ 10 26 WHz -1 . Sources with low 
apparent radio luminosity and slow speeds are unlikely to be powerful highly 



1 http://www. physics. purdue.edu/MOJAVE 
2 http://www. bu.edu/blazars/VLBAproject. html 
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relativistic jets oriented near the plane of the sky; nor can they be due to low 
luminosity high speed jets pointed close to the line of sight as the probability of 
these configurations is too low. 

With the extensive data base now available, many jets show clear curva- 
ture. A major question is to find out if the jet flow follows a preset curved 
trajectory or is ballistic. Ballistic motion can produce a curved jet if it origi- 
nates from a rotating nozzle, as might be induced if the AGN is part of a binary 
system. However, there is no clear model which explains all observations. In 
some sources like 0738 + 313, successive compo nents appear to follow the same 
curved trajectory as previo usly ejected f eatures ( Listerll2006l ) . However, in other 
sources, such as 1308+328, iListerl ( 20061 ) has shown that the observed pattern is 
consistent with that expected from a precessing nozzle having a 6.3 year period. 

Although gamma-ray emission is associated with flat spectrum radio blazars, 
the limited sensitivity and poor time resolution of the existing EGRET obser- 
vations has precluded a cl ear understanding of the relation between the radio 
and gamma-ray emission. iKovalev et all (|2005l ) found that gam ma-ray blazars 
have m ore compact structure than n on gamma-ray blazars, while Jorstad et al.l 
( 200ll ) and Kellermann et al. ( 2004) reported s o mewh at faster observed speeds 
for gamma-ray detected sources. Jorstad et al.l ( 200ll ) suggest that gamma-ray 
flares occur near the time of ejection of new super luminal components. 



4. Some sources of special interest 



M87 is one of the closest known relativistic jets and is of special interest since 
1 mas corresponds to a linear scale of only 0.08 pc and an angular v elocity of 1 



mas/y r corresponds to a linear velocity of 0.25 c. As early as 1964, IShklovskii 



(| 19641 ) argued that the jet in M87 was intrinsically two-sided but appeared one- 



sided due to differential Doppler boosting of the approaching and receding rel- 
ativistic jets. This was the first suggestion that bulk relativistic motion might 
play an i mportant role i n und erstanding extragalactic radio sources. However, 
curiously IKovalev et al. ( 20071 ) find no evidence at 2 cm over a period of 12 years 
for any motions in any of seven apparent fea tures greater t han a few percent c 
within about 1.5 pc of the jet base. However, iBiretta et al.l ( 19991 ) report super- 
luminal motion in t he optical feature kn own as HST-1 located 0.86 arcsec (60 pc) 
downstream, while ICheung et al.l (120071') report v/c ~ 4.5 i n the 20 cm radio fea- 
ture associated with HST-1. Also, Aharonian et al. ( 20061 ) have reported strong 
and variable TeV emission from M87. and ICheung et all (|2007ft argue that the 
TeV emission may originate in HST-1 and not at the base of the jet. Closer 
in to the base, 7 mm observations by iLv et al] ( 20071 ) with submilliarcsecond 

resolution suggests speeds which are a significa nt fraction of c . 

The high linear resolution 2 cm image of Kovalev et al. ( 20071 ) shows the 
beginning of the counter-jet. Kovelev et al. suggested that the apparent limb 
brightening of the main jet is due to a two component spine-sheath confi guration 
as expected from theoretical models (jAgudo et al.l l200ll ; iHarded 120071 ) . In the 
context of these models, the 2 cm observations see only the slowly moving outer 
sheath, while the faster m oving inner spine is responsib le for the observed high 
energy emission. However, iTavecchio Sz Ghisellinil ( 20081 ) speculate that the TeV 
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emission might be produced in the slowly moving outer layer surrounding a fast 
misaligned spine. 

3C 279 is one of the strongest ra dio blazers an d its j et kinematics have 



been studied for more than 30 years. iHoman et al.l (|2003l ) have reported an 



abrupt change in the trajectory of the brightest jet feature with an apparent 
chan ge in direction of ~ 26 deg and corresponding velocity change from 8c to 
13c. IHoman et al. (2003) concluded that the intrinsic bend was only about 1 
degree, but occured more than 1 kpc away from the base of the jet, and that 
the focusing and collimation of relativistic jets can occur well downstream. 

3C 111 is one of only a few known classi cal double lobed ra dio galaxies 
with observed superluminal motion in the core. iKadler et all (|2008l ) have stud- 
ied the evolution of the one sided radio jet in 3C 111, over a period of about 
10 years during a major flare and jet-plasma ejection, and found a complex 
pattern of separating components possibly indicative of forward- and backward 
directed shocks as well as trailing components formed in the wake of the primary 
perturbation. Such effects may be interpreted as a result of coupling to Kelvin- 
Helmholtz instability pinching modes from the interaction of the jet with the 
external medium. They complicate the interpretation of multi-epoch VLBI jet 
observations and provide a challenge for the descriptive analysis of jet structure 
as well as for numerical jet models." 



5. Summary and New Opportunities 

Although the parent population of AGN jets is mostly only mildly relativistic, in 
flux density-limited samples, Doppler boosting causes the observed distribution 
to be highly biased toward the fastest jets with observed Lorentz factors up 
to 35 or 40 which are aligned close to the line-of-sight. However, subluminal 
and stationary features are also observed. Although each jet appears to have a 
characteristic velocity which may reflect the bulk flow, quantitative analysis may 
be complicated by differences between the bulk plasma flow and observed pattern 
motions, by the apparent splitting and merging of individual features, by finite 
opening angles, and by transverse velocity gradients. On average quasars have 
larger Lorentz factors (7 ~ 8) than low luminosity AGN, and EGRET detected 
blazars may have somewhat larger Lorentz factors. The observed distribution 
of Lorentz factor and luminosity is consistent with Doppler boosting models, 
but there appear to be intrinsic differences as well as projection effects. Both 
ballistic motion characteristic of a rotating nozzle and curved trajectories with 
accelerations and decelerations are observed, and there is some evidence that in- 
situ acceleration along the jet may be important. The ejection of new features 
may be associated with radio or gamma-ray outbursts, but the evidence for this 
is inconclusive. Brightness temperatures up to 10 K are observed directly and 
are also inferred from the time scale of flux density variability. 

Over the course of the next few years, the bandwidth of the VLBA will 
be increased giving an improvement in sensitivity of up to a factor of 5 for 
sustained observations. In Europe, the EVN, especially when combined with 
eMERLIN, will give improved surface brightness sensitivity, and eVLBI will of- 
fer new target-of-opportunity possibilities, with real-time bandwidths of 1 Gbps. 
Continued radio observations, combined with the unprecedented GLAST all sky 



6 



Kellermann et al. 



monitoring, may give insight into the location and mechanism of gamma-ray 
outbursts. In 2012, VSOP-2 will give an order of magnitude improvement in 
both sensitivity and resolution over HALCA along with polarization capability. 
This will allow the transverse structures and magnetic field orientation in the 
M87 and other jets, to be studied in much better detail than previously pos- 
sible. In particular, VSOP-2 will better define the structure and evolution of 
jet-deflection/interaction surfaces and the role Kelvin-Helmholtz instabilities in 
sources such as 3C 111 and 3C 279. However, frequent time sampling and good 
u,v coverage will be needed to fully exploit these capabilities. 
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